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The electronic structure of artificial Mn atom arrays on Ag(lll) is characterized in detail with 
scanning tunnelling spectroscopy and spectroscopic imaging at low temperature. We demonstrate 
the degree to which variations in geometry may be used to control spatial and spectral distributions 
of surface state electrons confined within the arrays, how these are influenced by atoms placed within 
the structure and how the ability to induce spectral features at specific energies may be exploited 
through lineshape analyses to deduce quasiparticle lifetimes near the Fermi level. Through extensive 
comparison of dl/dV maps and spectra we demonstrate the utility of a model based upon two- 
dimensional s-wave scatterers for describing and predicting the characteristics of specific resonators. 

PACS numbers: 



I. INTRODUCTION 

A central aim of nanoscale science is to controllably 
modify material properties on a nanometre scale. One of 
its greatest achievements has been manipulation of sur- 
face atomic structure with the scanning tunnelling micro- 
scope (STM) , and some of the most striking demon- 
strations of this capability the construction of complex 
arrays of adatoms by the maneuvering and positioning of 
individual adatoms [Q . Synthesized nanoscale structures 
such as these, as well as other surface nanostructures 
including islands || [| and step arrays || act as res- 
onators for mobile surface electrons, with characteristics 
that may be controlled by the adjustment of geometrical 
parameters. 

Here we present a detailed analysis of the electronic 
structure of a particular class of resonator, surface 
nanocavities constructed from Mn atoms on Ag(lll). Of 
particular interest is the extent to which a simple scat- 
tering model, which takes into account the atomic struc- 
ture of the resonator, can successfully describe its spa- 
tial and spectral electron distributions. Recent work || 
has shown that the specific characteristics of individual 
resonators impact upon the electronic structure of ad- 
sorbates positioned within them, and may open up the 
possibility of atomic scale control of magnetic or chemical 
properties. 

In addition, we consider the use of these artificial 



nanoscale structures for probing quasiparticle lifetimes 
close to the Fermi energy. The lifetimes of excitations 
enter the description of many important surface phenom- 
ena such as the dynamics of charge and energy trans- 
fer, and for excitations in surface bands reflect a com- 
plex combination of decay processes via both surface- 
localised and extended states, with the additional influ- 
ence of screening that varies rapidly as the electron den- 
sity decreases outside the surface. Successful prediction 
of lifetimes remains an on-going theoretical problem , 
whilst measurements must overcome considerable chal- 
lenges (reviewed in Ref. ||) so that only recently |]] 
has a consistent account of the lifetime of a prototype 
quasiparticle, a hole in a surface band, been constructed, 
based upon tunnelling spectroscopy at the surface state 
band edge on noble metal (111) faces JTofl . The spatial 
decay of electron interference patterns observed in STM 
images provides an alternative approach first indicated 
by Avouris jy], and has been used to measure lifetimes 
well removed from the Fermi energy Ep ]l^ |. However, 
lifetime measurements for states close in energy to Ep 
are still lacking but are of particular interest because 
the relative contributions of phononic and electronic de- 
cay channels are expected to change significantly in this 
range. We indicate how electron resonators may be used 
for this purpose. 



II. METHODS 
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A. Experimental Procedures 

We study nanocavities constructed using Mn adatoms 
on the Ag (111) surface in a custom-built ultra-high 
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FIG. 1: (a) 60 x 60 A 2 constant current topograph of Mn 
adatoms on Ag(lll). At V = 512 mV and R = 731 MQ the 
Mn atoms occasionally exhibit tip-induced mobility (arrow). 
Image scanned at 1.9 s/line. (b) I /V spectrum acquired on 
top of a single Mn atom. Abrupt current fluctuations occur at 
V = 1.5 V and R — 1 GQ. Spectrum recorded at 133 mV/sec, 
starting at V = — 2 V. 



vacuum (UHV) STM, operating at a temperature of 
T = 4.6 K |l3f. The Mn atoms were evaporated onto 
a cold Ag substrate that had been prepared by stan- 
dard procedures in UHV. Electrochemically etched W 
tips were prepared in UHV by sputtering and annealing. 
Differential conductivity (dl/dV) spectra were recorded 
under open feedback loop conditions using a lock-in am- 
plifier with a sinusoidal modulation added to the sample 
bias. Moreover, maps of dl/dV were recorded simulta- 
neously with constant current imaging without opening 
the feedback. 

Manganese atoms on Ag(lll) image as 1.6 A high and 
11 A wide protrusions (full width at half maximum) at 
a sample voltage V = 1 V (Fig. |l]a), similar to previ- 
ous observations of adatom species on noble metal sur- 
faces 14 . Most of the STM images presented here 
have been recorded at V — 1 V and a tunnelling current 
I = 0.01 nA, i.e. a tunnelling resistance R = 100 Gf2; 
these values are chosen because the tip-sample interac- 
tion is sufficiently strong for R below 700 — 750 Mfl that 
atoms are occasionally moved during scans, as demon- 
strated in Figure [l]a. In I(V) spectroscopy of single Mn 
atoms we observe abrupt current changes at V « 1.5 V 
and R « 1 Gfi, Figure |l]b, indicating the onset of induced 
atomic mobility. To achieve acceptable signal-to-noise ra- 
tios in tunnelling spectroscopy elevated modulation am- 
plitudes were required (5 mV pp ). 

We construct artificial structures via controlled reloca- 
tion of individual Mn atoms, using a sliding process pio- 
neered by Eigler [Q . A Mn adatom is moved by placing 
the tip over the atom at large tip-sample separation (typ- 
ically R — 100 Gf2), reducing this separation (V = 3 V, 
R = 500 Mfi) and then moving the tip slowly across the 
surface (a rate of ~ one lattice constant per second) to 
the desired location. The adatom follows the tip motion, 
and is left in position by retracting the tip. 



B. Modelling 

To model the experiments it is important to take into 
account the constant-current mode of operation that is 
used in obtaining the dl/dV maps. The tunnelling cur- 
rent is calculated as Ml 



eV 



J(V,r)« J Ps (r,E)T{E,eV,z)dE, 



(1) 



where p s (r,E) is the local density of states (LDOS) of 
the sample at the position r = (x, y) of the tip and where 
the transmission probability T is given by 



T(E, eV, z) = exp 



4m 



eV - 2E) 



with tip-sample distance z and workfunctions (fit = 
4.55 eV and (fi s = 4.74 eV of tip and sample respectively 
M. Simulating the feedback loop, z was adjusted for 
each lateral position (z = z(r)) to ensure constant cur- 
rent /, which in turn was numerically differentiated to 
obtain dl/dV. In using (Q) we assume an electronically 
featureless tip, which is supported by the absence of any 
tip dependence to be seen during the experiments. 

For p s we consider an array of two-dimensional s- 
wave scatterers with coordinates Rj, j = 0, 1, . . . N — 1. 
This is the model introduced by Heller et al (H]]. We 
determine p s from the Green function 

p s (r,E) = -(l/n)ImG(r,r;E) (3) 

which is given by multiple-scattering theory in terms of 
the free-electron Green function fllql 



G (r,r';E) 



vm 
2h 



2 H^(K\r-r'\) 



(4) 



where k = y 2m*(E — Eo)/h 2 and the t-matrices that 
characterise the point scatterers: 

G(r,r';E) - G (r,r';E) 

+ Y / Go(r,R 3 ;E)T^ k Go(R k ,r';E) (5) 

j,k 

Ti\E) = t6 jk + tJ2Go(Rj,Rv,E)T lk (E). (6) 

A number of parameters enter the calculations. We 
use i-matrices corresponding to perfectly absorbing scat- 
terers, as identified by Heller et al. jlTj and used by 
them to describe the dl/dV spectrum of a 60 Fe atom 
corral on Cu(lll). This boundary condition reflects the 
strong coupling to bulk states |1!| that occurs at the 
adatoms. The effective mass m* = 0.42 m e and surface 
state band edge E = —67 meV have been determined 
previously |3, ^0|. The positions of the atoms and the 
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FIG. 2: Array of 28 Mn atoms forming a rectangle of size w 90 x 100 A 2 . Panels (a)-(h) show constant current topographs 
recorded at the indicated sample voltages. Panels (i)-(p) are dl/dV maps recorded simultaneously. 



lateral position of the tip are crucial in determining the 
exact shape of dl/dV spectra. For example, only the 
highest symmetry eigenstates contribute at the centre of 
a circular structure, but slightly off centre or if the shape 
is not exactly regular a mixing of state occurs, giving ad- 
ditional structure. Uncertainty enters these coordinates 
through the piezo calibration and possible anisotropy of 
the piezo sensitivity in the fast (x) and slow (y) scan- 
ning directions. Therefore an optimisation scheme was 
used based upon a % 2 fit using the Levenberg-Marquardt 
method in which the relative coordinates were fixed but 
independently scaled in the x and y directions. After 
the initial determination, the piezo calibration parame- 
ters were held fixed in all subsequent calculations. In a 
second stage, the calculations have been extended further 
to include a n ener gy dependent lifetime of the electronic 
states (Sect. HID). This was accomplished via a param- 
eterised imaginary self energy iE(E) added to the energy 
E in computing G from equations The impor- 

tance of the self energy in determining the peak ampli- 
tudes and widths in the LDOS has been noted previously 
|jl9| , Here we exploit this dependence to extract E 
from experimental data. 



III. RESULTS 

A. 28 atom rectangle 

Figure H displays data from an array of 28 Mn adatoms. 
Images a-h are constant current topographs covering a 
range of different voltages, recorded at values indicated in 
the figure. The adatoms, which appear as bright protru- 
sions in every case, have been arranged into a rectangular 
geometry with approximate dimensions of 90 x 100 A 2 . 
Inside this structure a weak corrugation is found, orig- 
inating from the confinement of the surface state elec- 
trons. A nonlinear height scale has been used for dis- 
play, in order to render this weak corrugation easily dis- 
cernible; this has the side effect that the adatoms appear 
rather blunt. At negative voltages (i.e. when tunnelling 
from occupied levels in the sample) the constant current 



topographs exhibit a rectangular pattern of four maxima 
within the confining array of Mn atoms, which is still 
discernible at voltages well below the bottom of the sur- 
face state band at Eq — —67 meV. This pattern persists 
owing to the preferential tunnelling from occupied states 
at the Fermi level. At positive voltages, imaging unoccu- 
pied states, the pattern becomes more complex and varies 
with V, the most significant tunnelling now occuring at 
energies corresponding to the varying bias voltage. 

Images i-p show dl/dV maps that were recorded si- 
multaneously with the topographic scans. Below Eq, 
the adatoms cause clear features, whereas no corruga- 
tion is discernible inside the array. In contrast, above Eq 
a clear wave pattern arises within the structure, while the 
adatoms are only weak features which conveniently serve 
to indicate the rectangular boundary. The evolution with 
energy of the internal pattern is as expected, represent- 
ing a progression through the nodal structures of the 
lowest lying "particle-in-a-box" eigenstates of the rectan- 
gular cavity, \i/) n ,m(x, y)\ oc | sin(nirx / X) sin(m7r?//y)| 2 
where X and Y are the linear dimensions of the struc- 
ture. At low energy, Fig. a single maximum oc- 
curs (n, m = 1, 1). Next, two maxima appear along the 
longer axis of the rectangle (Fig. |^k), corresponding to 
the n, m = 1,2 level. The —10 mV map (Fig. ||l) show- 
ing four maxima, two along each axis of the rectangle, 
reflects contributions from both the second (n,m = 1,2) 
and third (n,m = 2,1) levels. The states are in real- 
ity resonances that in general overlap in energy, so that 
dl / dV maps can contain contributions from several. At 
higher energies, the pattern becomes increasingly com- 
plex. 

Calculated results are displayed in Fig. |^ where pan- 
els a-g are topographs, h-n are dl/dV maps and o-u are 
LDOS maps. Close inspection and comparison with Figs. 
P>h, and j^j-p reveals that the experimental features are 
reproduced well by the model calculations. As pointed 
out previously pQ ], experimental dl/dV data recorded 
by the methods used here are affected by the variations 
in tip-sample separation which give rise to the corruga- 
tion in the topographs. Consequently dl/dV maps are 
not directly comparable to LDOS maps, a fact clearly 
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visible by comparing Figs. ||l and s, or and t, al- 
though the similarity between the two tends to increase 
with increasing bias voltages. Indeed, around zero bias 
it is topographs and not dl/dV that are are expected to 
resemble the LDOS, and this is clearly the case in Fig. 

I 

A more detailed comparison between experiment and 
theory is possible in the case of dl/dV spectra. An exam- 
ple measured near the centre of the rectangular resonator 
is shown in the left panel of Fig. ^. Tunneling from a 
s-state tip directly at the centre of an ideal rectangu- 
lar structure can only take place into surface state levels 
with odd quantum numbers n, m. This is the origin of 
the large peaks at —50 meV and between and 40 meV. 
Irregularities in the structure and deviations in position 
of the tip from the the exact centre introduce additional 
weak tunneling channels, e.g. near —25 meV, whilst the 
intrinsic width of the confined levels further modifies the 
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FIG. 4: Left: Measured dl/dV spectrum and calculated 
LDOS of the 28 atom rectangle. Right: dl /dV spectrum and 
theoretical LDOS of a 220 A diameter ring of 35 Mn adatoms. 



spectrum. Comparing the measured spectrum with the 
LDOS calculated using the two-dimensional scattering 
model that includes the structural irregularities and cor- 
rect tip position, displayed in the same panel in Fig. |J, 
we see the latter faithfully reproduces the differential con- 
ductivity spectrum. Note that the characteristic temper- 
ature associated with Kondo behaviour in Mn/Ag alloys 
is Tk — 40 mK [^2) and thus does not affect the line 
shape p3[ under the present experimental conditions. 



B. 35 atom ring 

We have also performed similar experiments and calcu- 
lations on other geometries, such as adatom rings. The 
right panel of Figure |5] contains a comparison between 
the differential conductivity spectrum, measured approx- 
imately 5 A off centre in a 220 A diameter ring formed 
from 35 Mn adatoms, and the corresponding LDOS ob- 
tained using the s-wave scattering model. The spectrum 
contains a large number of peaks, the off-centre position 
enabling low symmetry eigenstates to provide channels 
for tunnelling, but these are reproduced in the theoreti- 
cal LDOS spectrum calculated using the precise atomic 
positions of the Mn ring. We conclude that the s-wave 
scattering model provides an adequate basis for mod- 
elling the spatial and spectral electron distributions in 
electron resonators. 



C. Adatoms in resonators 

Having characterised the properties of electron res- 
onators we have next investigated the interaction with 
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FIG. 5: Topographs of a 28 atom rectangle with a central 
Mn adatom recorded at the indicated sample voltages (a-c, 
top row), and dl/dV maps of the same structure with (d-f, 
middle row) and without (g-i, bottom row) the central atom. 

an additional atom introduced into the nanocavity re- 
gion. These structures were formed by maneuvering one 
of the boundary Mn atoms into the central region, and 
then moving a Mn atom elsewhere on the surface into the 
vacant position in the boundary, restoring the integrity 
of the resonator. 

In Figure || we show dl/dV maps of a rectangular 
structure made from 28 Mn atoms and measured at three 
different voltages, contrasting the images obtained with 
and without a Mn atom inside the structure. At first 
glance the presence of the central atom appears to barely 
affect the wave pattern inside the rectangle, apart from 
in the immediate vicinity of the centre atom. There, 
the difference is striking. At voltages of V — —10 mV 
and V = 50 mV where dl/dV is at a minimum in the 
empty structure (Figs. ||g and i), dl/dV exhibits a peak 
in the centre upon insertion of the atom (Figs. |^d and f). 
At V — 10 mV, however, the maximum observed in the 
empty rectangle (Fig. |B|h) changes to a minimum (Fig. 

This "contrast reversal" has been previously reported 
and explained as a consequence of the strong hybridi- 
sation between states on the Mn atom and substrate 
surface electron levels ||, and has been reproduced in 
calculations of the local density of states directly on 
Mn atoms in resonators of various geometries using full 
three-dimensional multiple scattering calculations. Here 
we consider how this interaction extends away from the 
central Mn atom and the degree to which the two- 



FIG. 6: Evolution of recorded dl /dV spectra with increasing 
lateral distance from the centre of a 16 Mn atom ring radius ~ 
40 A (a) without and (b) with a central Mn atom, and spectra 
calculated using the two-dimensional scattering model of the 
LDOS for the same systems, (c) without and (d) with the 
central atom. The vertical dashed lines mark the energies of 
the main spectral features present at the centre of the empty 
structure. 



dimensional scattering model describes the behaviour. 
Fig. U shows a series of dl/dV spectra measured at var- 
ious distances from the centre of an empty adatom ring, 
and when a Mn atom has been introduced. Also shown 
arc the LDOS calculations for identical geometries. Com- 
paring dl/dV spectra in Figs. ||a and b at central posi- 
tions closest to the location of the new adatom, it is clear 
that spectral peaks in the empty resonator coincide with 
troughs of the occupied structure (vertical dashed lines). 
This is evidence for strong hybridisation Q . Moving lat- 
erally away from this central position the characteristic 
peak structure of the empty resonator returns, with the 
exception of the lowest peak near —40 mV which remains 
absent even at large displacements. This prominent fea- 
ture is associated with a nodeless resonator eigenstate 
that has its maximum at the centre of the ring, and so is 
particularly sensitive to the introduction of a Mn atom 
into the central region. 

The two right-most panels in Fig. || show the LDOS 
calculated for this system using the two-dimensional scat- 
tering model, at lateral positions at which dl /dV spectra 
were recorded. The prominent features in the experimen- 
tal spectra are reproduced well in the theoretical curves, 
again showing that the simple scattering model, using 
the exact atomic locations of the Mn atoms making up 
the structure, gives an accurate and valid account of the 
electronic structure of resonators. The calculations also 
confirm that the dominant state in the empty resonator 
is absent in the spectrum of states of the occupied res- 
onator, so that the insensitivity of dl/dV maps to the 
presence of the central atom away from its immediate 
vicinity that was suggested by Fig. |B| clearly does not 
extend to all sample voltages. Note that in Fig. ||d we 
have omitted theoretical spectra within 10 A of the cen- 
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FIG. 7: Effect of a constant imaginary part E of the self 
energy on the LDOS of a 36 atom circle with 210 A diameter, 
(a) Calculated LDOS at the centre for various E indicated in 
the figure. Variation of peak (b) heights and (c) widths with 



tral Mn atom. Within this radius (which is much greater 
than the topographic radius of the Mn atoms) the dif- 
ferential conductivity also has large contributions from 
channels involving bulk states scattered by the adatom 
p4j , and these are not included in the model. Also omit- 
ted in the model is any spatial information about the 
adatom orbitals that would also be required to model 
dl/dV spectra (and maps) in its immediate vicinity. 



D. Energy dependent self energy 

The basic idea of our approach to lifetime measure- 
ments is straightforward [p.3| . Exploiting the properties 
of the electron resonators identified in the previous sec- 
tions we construct nanocavities with electronic features 
in the vicinity of Ep. The lineshape of the spectral fea- 
tures are affected by lifetime-limiting processes so that 
through analysis the lifetime r can be determined. By 
changing the geometry of nanoscale structures it is pos- 
sible to tune the eigenstate energies of the structures, 
and thus generate states at various energies in order to 
obtain the energy dependence of the lifetime. The two- 
dimensional scattering model validated above enables the 
design of suitable structures. We note that a fully quan- 
titative analysis will also require information on the en- 
ergy dependence of the tip density of states, which is not 
expected to be constant over a wide energy range, and 
likewise the bulk density of states of the crystal which 
causes a background signal in dl/dV spectra. Here, we 
seek to demonstrate the concept and so neglect these in- 
fluences. Following Li et al. |10| we describe the lifetime 
effects in terms of an effective imaginary part E of the 
self energy, where r = ?i/(2E). 

The effect of E on the LDOS of an artificial nanoscale 
resonator is shown in Fig. Qa for a 210 A diameter 
ring. The calculation for vanishing E, i.e. infinite life- 
time, yields a series of sharp resonances (Fig. 0a, top 
curve). Note that the spectrum corresponds to the centre 
of the structure. The high symmetry gives the "cleanest" 
spectrum, avoiding difficulties that exist at off-centre po- 
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FIG. 8: E(i5) extracted from analysis of dl /dV spectra of a 
115 A radius 35 atom circular resonator (diamonds, squares) 
and a 40 A radius 16 atom circular resonator. The E values 
are plotted at the position of the respective spectral features. 
The inset shows the equivalent lifetimes r = ft/(2E) of hot 
holes and electrons. 



sitions where resonances can overlap. With increasing E 
there is a reduction of the peak heights and a concomi- 
tant broadening, especially apparent for the lowest of the 
states which have the sharpest intrinsic profile. These 
trends are quantified by fitting a series of Lorentzians to 
the data, and Figure |7p and c display the widths and 
height extracted through this analysis. As E increases, 
the peak heights decrease monotonically, the lowest (in 
energy) peaks most rapidly, whilst there is an essentially 
linear dependence of the peak width on E. 

Figure || shows values of E(i?) obtained from an anal- 
ysis of measurements from two nanostructures. For each, 
X 2 fits were made over limited energy ranges (~ 100 
meV) adjusting coefficients of a linear parameterisation 
of E. We estimate the error of E thus determined to 
be of the order of ±30 %. The diamonds and squares 
are derived from dl/dV spectra for a 35 atom circu- 
lar resonator, and the filled circles were obtained from 
a 16 atom circle. Below the Fermi energy (E — 0) E ap- 
proaches a value obtained previously from spectra taken 
on clean terraces which yielded E = 3 meV g|. There 
is a clear increase of E at large energies, and we find 
E(_E = 0.5 eV) ~ 20 meV, which is comparable with an 
extrapolation down to this energy of data for the Ag(lll) 
surface state extracted from the spatial coherence of in- 
terference patterns at energies E > 1 eV by Biirgi et 
al. |fl2|] . Finally, at energies close to Ep we observe a 
marked decrease of E which corresponds to an increased 
lifetime of the corresponding quasiparticles (inset in Fig. 
|8|). Increased lifetimes are indeed expected as both the 
phononic and electronic decay channels become less ef- 
ficient. While further data close to Ep is necessary to 
quantitatively characterize this effect the first results are 
encouraging. Future work will also address the relation 
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of the self-energies determined from atomic arrays to that 
of an unperturbed surface. We note that data from res- 
onators of different diametres (40 A and 115 A) yield 
similar S(-E). This similarity appears to indicate that 
the influence of the array may be rather small. 

IV. SUMMARY 

Using scanning tunnelling spectroscopy and spectro- 
scopic imaging at low temperature we have investigated 
artificial Mn atom arrays on Ag(lll). These act as elec- 
tron resonators and we have demonstrated the degree 
to which variations in geometry may be used to control 
spatial and spectral distributions of confined surface state 



electrons. This control may be exploited in manipulating 
the local density of states of adsorbates introduced into 
the resonator, or as we have demonstrated, to generate 
spectral features at specific energies in order to perform 
lineshape analyses from which quasiparticlc lifetimes may 
be deduced. Extensive comparison of dl/dV maps and 
spectra indicated the utility of a model based upon two- 
dimensional s-wave scatterers for describing and predict- 
ing the characteristics of specific resonators. 
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